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ABSTRACT 



(N 
> 

, In classical gravity deviations from the predictions of the Einstein theory are often discussed within 

the framework of the conformal Newtonian gauge, where scalar perturbations are described by 
I two potentials (p "0- In this paper we use the above gauge to explore possible cosmological 

consequences of a running Newton's constant as suggested by the nontrivial ultraviolet 

fixed point scenario arising from the quantum field-theoretic treatment of Einstein gravity with 
a cosmological constant term. Here we focus on the effects of a scale-dependent coupling on the 
so-called gravitational slip functions rj = — 1, whose classical general relativity value is zero. 
Starting from a set of manifestly covariant but nonlocal effective field equations derived earlier, 
we compute the leading corrections in the potentials (p and for a nonrelativistic, pressureless 
fluid. After providing an estimate for the quantity ??, we then focus on a comparison with results 
obtained in a previous paper on matter density perturbations in the synchronous gauge, which 
gave an estimate for the growth index parameter 7, also in the presence of a running G. Our 
results indicate that, in the present framework and for a given G(n), the corrections tend to be 
significantly larger in magnitude for the perturbation growth exponents than for the conformal 
Newtonian gauge slip function. 
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1 Introduction 



Recent years have seen the development of a fascinating variety of alternative theories of gravity, 
in addition to the more traditional alternate frameworks, which used to include just Brans-Dicke, 
higher derivative, effective quantum gravity and supergravity theories. Some of the new additions 
to the by now rather long list include dilaton gravity, f{R) gravity, torsion gravity, loop quantum 
gravity, holographic modified gravity, and a few others, just to cite here a few representative 
examples. All of these theories eventually predict some level of deviation from classical gravity, at 
short- or long-distance scales, which is often parametrized either by a suitable set of post-Newtonian 
parameters, or more recently, by the introduction of a gravitational slip function [HE]- 

In this paper, we will focus on the analysis of departures from general relativity (GR) in the 
gravitational slip function, obtained in the framework of the conformal Newtonian gauge, and within 
the rather narrow context of the nontrivial ultraviolet fixed point scenario for Einstein gravity with 
a cosmological term. Thus, instead of looking at deviations from GR at very short distances, due 
to new interactions such as the ones suggested by string theories p], we will be considering here 
infrared effects, which could manifest themselves at very large distances. 

The specific nature of the scenario we will be investigating here is motivated by the field- 
theoretic treatment of models for quantum gravity, based on the (minimal) Einstein action with 
a bare cosmological term. The theory's long-distance scaling properties used as the basis for the 
present work follow from the existence of a nontrivial ultraviolet fixed point of the renormalization 
group in Newton's constant G. The latter is inaccessible by direct perturbation theory in four 
dimensions, and can be shown to radically alter the short- and long-distance behavior of the theory 
when compared to more naive, perturbative expectations. The renormalization- group origin of 
such fixed points was first discussed in detail by Wilson some time ago for scalar and self-coupled 
fermion theories [4J . The general field-theoretic methods were later extended and applied to gravity, 
where they are now referred to as the nontrivial UV fixed point, or asymptotic safety, scenario 
It is fair to say that so far this is the only field-theoretic approach known to work consistently 
in other not perturbatively renormalizable theories, such as the nonlinear sigma model above two 
dimensions [B]. While perhaps still a bit mundane in the context of gravity, such nontrivial fixed 
points are well studied and well understood in statistical field theory, where they generally describe 
phase transitions between ordered and disordered ground states, or between weakly coupled and 
condensed states. 
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The paper is organized as follows. First (Sec. 2) we recall the effective covariant field equations 
describing the running of G, and describe briefly the nature of various objects and parameters 
entering the quantum nonlocal corrections; a more complete description of the basic setup can 
be found in our previous papers on the subject, and will not be repeated here. We then discuss 
the zeroth order (in the metric fluctuations) field equations and energy-momentum conservation 
equations for the standard homogeneous isotropic metric, with a running G(n). Later (Sec. 3) 
we extend the formalism to deal with small metric and matter perturbations, and list the relevant 
field and energy conservation equations to first order in the perturbations in the comoving gauge. 
These above results are then (Sec. 4) reexpressed in two other choices of gauge, the synchronous 
and the conformal Newtonian gauge. The latter choice of gauge allows us to extract an expression 
for the gravitational slip function rj due to G(n) (Sec. 5). This quantity is then evaluated within 
the context of a AC DM model, for redshifts corresponding to the present era (z=0). The resulting 
correction is then compared to current astrophysical observations, as well as to our previous results 
(and observations) regarding the corrections due to (?(□) to the matter density perturbation growth 
exponents. The conclusions provide an interpretation of the theoretical results, and their associated 
uncertainties, in view of present and future high precision determination of the gravitational slip 
function and growth exponents. 



2 Running Newton's Constant G(n) 



As mentioned in the introduction, it is not the purpose of this paper to provide a satisfactory 
description, or motivation, for the running of G that arises in the quantum-field-theoretic treatment 
of Einstein's gravity with a cosmological term. Here we only provide a brief summary, and only 
the most relevant formulas will be given for later reference; a more complete set of references can 
be found, for example, in [7J. 

The running of Newton's constant G has been computed both on the lattice in four dimensions 
[51 [S] , and in the continuum within the framework of the background field expansion applied to 
d = 2+e spacetime dimensions [5l[T0], and later also using truncation methods applied in d = 4 [11 j. 
In either case one obtains a momentum-dependent G{k'^), which eventually needs to be reexpressed 
in a suitable coordinate-independent way, so that it can be consistently applied to more general 
problems, involving arbitrary background geometries. The first step in analyzing the consequences 
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of a running of G is therefore to rewrite the expression for G{k'^) in a coordinate-independent way, 
either by the use of a nonlocal Vilkovisky-type effective gravity action |12l [T3] , or by the use of a 
set of consistent effective field equations. In going from momentum to position space one employs 
k"^ — )• — □ , which then gives for the quantum-mechanical running of the gravitational coupling the 
replacement G — t- G{0). Then the running of G is given in the vicinity of the UV fixed point by 



G{a) = Go 



1 \l/2i^ 

l + co(^j +... 



(2.1) 



where □ = g^'^'V^Vi, is the covariant d'Alembertian, and the dots represent higher order terms 
in an expansion in l/(^^n). Current evidence from Euclidean lattice quantum gravity points 
toward cq > (implying infrared growth) and ~ ^ [9j. Within the quantum- field-theoretic 
renormalization-group treatment, the quantity ^ arises as an integration constant of the Callan- 
Symanzik renormalization-group equations. 

One issue of great relevance to the physical interpretation of the results, is therefore a correct 
identification of the renormalization-group invariant scale ^. A number of arguments, mostly based 
on nonperturbative lattice results and scaling considerations involving the gravitational Wilson loop 
and its relevance for large scale observable curvature can be given in support of the suggestion 
that the dynamically generated infrared cutoff scale ^ (analogous to the A-^g^ of QCD) can be quite 
large in the case of gravity (for a recent review, see Ref. [7|). These arguments would then suggest 
that the new scale is naturally expected to be related to the large scale average curvature, and 
thus could be of cosmological relevance. 



A A . (2.2) 



These considerations then lead to a more concrete quantitative estimate for the scale in the running 
G{0) of Eq. (j2.ip . namely ^ ^ ~ 1.51 x 10^'^cm. Moreover, from these types of arguments 



one would also infer that the constant Gq in Eq. p.ip can, to a very close approximation, be 
identified with the laboratory value of Newton's constant, -v/Gq ~ 1.6 x 10^'^'^ cm. The running of 
G envisioned above would then remain in agreement with laboratory and solar system precision 
tests of general relativity. 

The appearance of the d'Alembertian □ in the running of G naturally leads to both a nonlocal 
effective gravitational action and a corresponding set of nonlocal modified field equations. In the 
simplest scenario, instead of the ordinary Einstein field equations with constant G 

Riiu - \g^iuR + >^gtJiv = ^-^GT^y , (2.3) 
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one is now led to consider the modified effective field equations 

R^^u - \ gi,vR + Xg„u = 8^G(n)r^^ (2.4) 

with the nonlocal term due to the G{0). By being manifestly covariant, these equations still satisfy 
some of the basic requirements for a set of consistent field equations incorporating the running of G. 
Not unexpectedly though, the new nonlocal equations are much harder to solve than the original 
classical field equations for constant G. 

The effective nonlocal field equations of Eq. (j2.4p can be recast in a form very similar to the 
classical field equations, but with a new source term T^i, = [G(n)/Go] Tf^w defined as the effective, 
or gravitationally dressed, energy- momentum tensor 116) . Ultimately the consistency of the 
effective field equations demands that it be exactly conserved, in consideration of the contracted 
Bianchi identity satisfied by the Ricci tensor. In this picture, therefore, the running of G can 
be viewed as contributing to a sort of vacuum fluid, introduced in order to account for the new 
gravitational quantum vacuum-polarization contribution. 

Due the appearance of a negative fractional exponent in Eq. (|2.ip . the covariant operator ap- 
pearing in the expression for G{0) has to be suitably defined by analytic continuation. This can 
be done, for example, by computing for positive integer n, and then analytically continuing to 
n — )• —l/2iy |15] . Equivalently, G(n) can be defined via a regulated parametric integral represen- 
tation [17], such as 

/ 1 \ l/^'^ 1 /-oo 

( A ? ) = -TT-: / da aV2.-i ^-a (2.5) 

V-D(5) + ^2; n^) Jo 

where ^ — )■ is a suitable infrared regulator. As far as the calculations in this paper are concerned, 
it will not be necessary to commit oneself to an unduly specific form for the running of G{0). 
Thus, for example, although the lattice gravity results only allow for a nondegenerate phase for 
the case cq > 0, it will nevertheless be possible later to have either sign for the correction in 
Eq. (j2.1|) . We note here that a running cosmological constant X{k) — t- A(n) causes a number of 
mathematical inconsistencies [18] within the manifestly covariant framework, described here by 
the effective field equations of Eq. ()2.4p . Indeed if one assumes for the running part of A(n) ~ 
(^^□)~'^, then the infrared regulated expression in Eq. (|2.5p gives no running of A, after using 
the identity V \g^u = 0. This last conclusion is in general agreement with the field-theoretic 



^ To be a bit specific, consider the case of a scale dependent A(fc), wliicli we will write here as A = Ao + 5\{k). 
Let us also assume, for concreteness, that 5\{k) ^ ci(fc'^)~'^, where ci and a are some constants, and then make the 
transition to coordinate space by using P -a. Thus <5A(D) ~ (-□ + ^i'^y, where one should be careful and use 
the infrared regulated expression in Eq. I|2.5p . The effective field equations will then contain a term ^ SX{0) ■ 

= I ci Jjj°° da ci'~^ e~" ("'-'(9)+'^ ) . g^^ — 1 ci {^J?)'" ■ Qiiv, which is still gives just a constant multiplying the 
metric g^iv 
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results of the nontrivial renormalization-group fixed point scenario [7j, thereby providing perhaps 
an independent consistency check. Note that this rather general argument also applies to possible 
additional contributions from nonzero vacuum expectation values of matter fields, such as the 
Higgs. As a result, in the present quantum-field-theoretic motivated framework A is assumed not 
to run. 



2.1 Zeroth order effective field equations with G(n) 

A scale-dependent Newton's constant is expected to lead to small modifications of the standard 
cosmological solutions to the Einstein field equations. Here we will summarize what modifica- 
tions are expected from the effective field equations on the basis of The starting point is 
the quantum effective nonlocal field equations of Eq. (j2.4p . with defined in Eq. (|2.ip . In 
the Friedmann-Lemaitre-Robertson- Walker (FLRW) framework these are applied to the standard 
homogeneous isotropic metric 

dT^ = dt^ - a^{t) I ^ '^'^^^^ + (de'^ + sin^ 6 dip"^^ | k = 0,±l. (2.6) 

In the following, we will only consider the case k = (spatially flat universe). It should be noted 
that there are in fact two related quantum contributions to the effective covariant field equations. 
The first one arises because of the presence of a nonvanishing cosmological constant A ~ 3/^^, 
caused by the nonperturbative quantum vacuum condensate < ii > 7^ [14J. As in the case of the 
standard FLRW cosmology, this is expected to be the dominant contributions at large times t, and 
gives an exponential (for A > 0), or cyclic (for A < 0) expansion of the scale factor. The second 
contribution arises because of the explicit running of G(n) in the effective field equations. 

The next step therefore is a systematic examination of the nature of the solutions to the full 
effective field equations, with G{0) involving the relevant covariant d'Alembertian operator 

□ = 5^"V^V, (2.7) 

acting on second rank tensors as in the case of T^u- To start the process, we will assume that T^jy 
is described by the perfect fluid form, 

T^u = [pit) + p{t) ]u^,uy + gf,^ p{t) (2.8) 

for which one needs to compute the action of on T^^,, and then analytically continues the 
answer to negative fractional values of n = —l/2u. The results of |15l 116^ ITT l I18j then show that a 
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nonvanishing pressure contribution is generated in the effective field equations, even if one initially 
assumes a pressureless fluid, p{t) = 0. After a somewhat lengthy derivation one obtains for a 
universe filled with nonrelativistic matter {p=0) the following set of effective Friedmann equations. 



+ 



a^t) ' a2(t) 

for the tt field equation, and 
k 



8TTG(t) A 



1 



ctit/0 



l/u 



2d(t) 



Stt Gf) 



ct{t/0 



l/u 



+ 



P(0 + 3 



p[t) + A 



(2.9) 



(2.10) 



a\t) a\t) ' ait) 3 

for the rr field equation. In the above expressions, the running of G appropriate for the Robertson- 
Walker metric is 



1 + ct 



+ ... 



(2.111 



with ct of the same order as cq in Eq. (I5.48P [15j (in the quoted reference the estimate q ~ 0.450 cq 
was given for the tensor box operator). From the above form of 6G{t) one sees that the amplitude 
of the quantum correction is actually proportional to the combination cq/^^ for u = 1/3. Note also 
that the running of G induces an effective pressure term in the second (rr) equation, due to the 
presence of a relativistic fluid whose origin is in the vacuum-polarization contribution. Another 
noteworthy general feature of the new field equations is the additional power-law acceleration 
contribution, on top of the standard one due to the A term. 



2.2 Introduction of the Wvac parameter 



It was noted in [15\ 118] that the field equations with a running G, Eqs. (|2.9p and (j2.10p . can 
be recast in an equivalent, but slightly more appealing, form by defining a vacuum-polarization 
pressure p^ac and density pvac, such that for the FLRW background one has 

. ^ 6G(t) , , . ^ 1 SG(t) , , , , 

Pvacit) = -p^ Pit) p.acit) = - Pit) . (2.12) 

From this viewpoint, the inclusion of a vacuum-polarization contribution in the FLRW framework 
seems to amount to a replacement 

pit) ^ pit) + pvacit) Pit) ^ Pit) + Pvacit) (2.13) 
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in tlie original field equations. Just as one introduces the parameter w, describing the matter 
equation of state, 

p{t) = wp{t) (2.14) 
with w = for nonrelativistic matter, one can do the same for the remaining contribution by setting 

it) . (2.15) 

We should remark here that the original calculations |15j . and more recently [T^ which included 
metric perturbations, also indicate that 

Wyac = I (2.16) 

is obtained generally for the given class of G{0) considered, and is not tied therefore to a specific 
choice of u, such as u = |. 

The previous, slightly more compact, notation allows one to rewrite the field equations for the 
FLRW background in an equivalent form, which we will describe next. We note here that, when 
dealing with density perturbations, we will have to distinguish the background, which will involve a 
background pressure (p) and background density (p), from the corresponding perturbations which 
will be denoted here by 6p and 5 p. With this notation and for constant Go, the FLRW field 
equations for the background are written as 

3^ = 8^Gop(t) + A 

^ + 2^ = -SttGoMO + A. (2.17) 
a^{t) a[t) 

Then in the presence of a running G(n), and in accordance with the results of Eqs. ()2.9p and ()2.10p . 
the modified FLRW equations for the background read 

\2 



a^(t) a[t) V Go / 

using the definitions in Eqs. ()2.14p and (12.150 . here with pvacif) = WvacPvac{t). 

Of course the procedure of defining a p^ac ^-nd a p^jac contribution, arising from quantum vacuum- 
polarization effects, is not necessarily restricted to the FLRW background metric case. In general 
one can decompose the full source term in the effective nonlocal field equations of Eq. (|2.4p . making 
use of 



as two contributions, 

^ G(a) T,, = {l + T,, = T,, + T;T . (2.20) 
The latter involves the nonlocal part 

= (2.21) 

Consistency of the full nonlocal field equations requires that the sum be conserved, 

V'^(V+T;^'=) =0. (2.22) 

In general one cannot expect that the contribution T™'^ will always be expressible in the perfect 
fluid form of Eq. (j2.8p . even if the original T^y for matter (or radiation) has such a form. The 
former will in general contain, for example, nonvanishing shear stress contributions, even if they 
were originally absent in the matter part. 



3 Relativistic treatment of matter density perturbations 

Besides the modified cosmic scale factor evolution just discussed, the running of G(n), as given 
in Eq. (|2.ip . also affects the nature of matter density perturbations on large scales. In computing 
these effects, it is customary to introduce a perturbed metric of the form 

dr^ = dt^ - G? {bij + hij) dx'dx^ , (3.1) 

with a{t) the unperturbed scale factor and hij{x,t) a small metric perturbation, and /iqo = /lio = 
by choice of coordinates. After decomposing the matter fields into background and fluctuation 
contribution, p = p + 6p, p = p + 6p, and v = v + 6v, it is customary in these treatments to expand 
the density, pressure and metric perturbations in spatial Fourier modes, 

6pi^,t) = 6p^it)e'^-^ 6p{^,t)=6p^{t)e'^-^ 

<^v(x,t) = 5vq(t)e^'i-" hij{^,t) = h^i,{t)e'^-^ (3.2) 

with q the comoving wave number. Once the Fourier coefficients have been determined, the original 
perturbations can later be obtained from 

Mx,t) = / j0j^e-^''"'SPnit) (3.3) 
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and similarly for the other fluctuation components. Then the field equations with a constant Gq 
[Eq. (j2.3p ] are given to zeroth order in the perturbations by Eq. ()2.17p . which fixes the three 
background fields a{t), p{t), and p{t). Note that in a comoving frame the four- velocity appearing 
in Eq. ()2.8p has components = 1, = 0. Without G'(n), to first order in the perturbations and 
in the limit q — t- the field equations give 

ait) 



ait) 



hit) = Stt Go pit) 6it) 



hit) + 3^ hit) = - 247r Go u; (3.4) 
ait) 

with the matter density contrast defined as = 5 pit) /pit), hit) = huit) the trace part of hij, 
and 1/7 = for nonrelativistic matter. When combined together, these last two equations then yield 
a single equation for the trace of the metric perturbation, 

hit) + 2^ hit) = -8^Go(l + 3/«)p(t)5(t) . (3.5) 
ait) 

From first order energy conservation, one has — ^ (1 + w) hit) = Sit), which then allows one to 
eliminate hit) in favor of 5(t), which then allows one to obtain a single second order equation for 
the density contrast (5(t). In the case of a running G(n), the above equations need to be rederived 
from the effective covariant field equations of Eq. (j2.4p . and lead to several additional terms not 
present at the classical level )18] . 



3.1 Zeroth order energy-momentum conservation 



As a first step in computing the effects of density matter perturbations, one needs to examine 
the consequences of energy and momentum conservation, to zeroth and first order in the relevant 
perturbations. If one takes the covariant divergence of the field equations in Eq. (j2.4p . the left- 
hand side has to vanish identically because of the Bianchi identity. The right-hand side then gives 
VA* (t^^ + T^^c^ = 0, where the fields in T^^^ can be expressed, at least to lowest order, in terms 
of the Pvac and p^ac fields defined in Eqs. (j2.12p and (j2.15p . The first equation one obtains is the 
zeroth (in the fluctuations) order energy conservation in the presence of G(n), which reads 



^ait) 
ait) 



^ M ^ ^Git) 

(1 + «;) + (1 + W^ac) -pT^ 



,W + ^,W+(l + ^|!))^(„.„. (3.6) 



In the absence of a running G these equations reduce to the ordinary mass conservation equation 
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for w = 0, 

It is often convenient to solve the energy conservation equation not for p{t), but instead for p(a). 
This requires that, instead of using the expression for G{t) in Eq. (j2.1ip . one uses the equivalent 
expression for G{a) 

G{a) = Go(l + , (3.8) 



Go 

which is easily obtained once the relationship between t and a{t) is known (see discussion later 
Note for example that the solution to Eq. (|3.6p can be written as 



, , da / 5G(a) SG'(a)\ , 
p{a) = const, exp <{ - / — ( 3 + + a ] \ . (3. 



3.2 Effective energy-momentum tensor involving pvac and pvac 

The next step consists in obtaining the equations which govern the effects of small field per- 
turbations. These equations will involve, apart from the metric perturbation hij, the matter and 
vacuum-polarization contributions. The latter arise from [see Eq. ()2.20p ] 

(l + T,, = T,. + (3.10) 

with a nonlocal TJ^^'^ = {SG{0)/Gq) T^^. Fortunately to zeroth order in the fluctuations the results 
of Ref. {15] indicated that the modifications from the nonlocal vacuum-polarization term could 
simply be accounted for by the substitution 

p{t) p{t)+p^ac{t) Pit) Pit) +pvacit) ■ (3.11) 

Here we will apply this last result to the small field fluctuations as well, and set 

SPdit) ^ Spciit) + dpc^vacit) SPciit) ^ Sp^it) + 6pqyacit) . (3.12) 

The underlying assumption is of course that the equation of state for the vacuum fluid still remains 
roughly correct when a small perturbation is added. Furthermore, just like we had pit) = w pit) 
[Eq. (j2.14p ] and pvadt) = w^ac Pvacit) [Eq. (j2.15p ] with w^ac = |, we now write for the fluctuations 

Spciit) =wSpc^it) 6pc^vacit) = Wyac^Pcivacit) , (3.13) 
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at least to leading order in the long wavelength limit, q — )• 0. In this limit we then have simply 

Sp{t) = W 6p{t) Spyac{t) = Wvac ^Pvacit) = Wyac ^^Jf^ ^P{t) , (3.14) 

with G{t) given in Eq. (j2.1ip . and we have used Eq. (|2.12p . now applied to the fluctuation 5p^ac{t), 

Spvacit) = ^-^6p{t) + ... (3.15) 

where the dots indicate possible additional 0{h) contributions. A bit of thought reveals that the 
above treatment is incomplete, since G{0) in the effective field equation of Eq. ()2.4p contains, for 
the perturbed Robertson- Walker metric of Eq. (|3.ip . terms of order hij, which need to be accounted 
for in the effective T^ac- Consequently the covariant d'Alembertian operator □ = g'^'^ V^V,/ acting 
here on second rank tensors, such as T^j^, 

(VuTa/s) = dfj_I^al3 — ^tfi^XalS — ^a/i^uXlS — ^'^fi^uaX ; (3.16) 

needs to be Taylor expanded in the small field perturbation hij, 

a{g) = + □(!)(/!) + 0(/i2) . (3.17) 



One then obtains for G(n) itself 
G{a) = Go 



1+^ i ) + 

^ ^i/u yaio) + □{i)(/i) + o(/j2) ' ^ • • • 



(3.18) 



which requires the use of the binomial expansion for the operator {A+By^ = A^'^ - A^'^ B A^'^ + . . .. 
Thus for sufficiently small perturbations it should be adequate to expand G(n) entering the effective 
field equations in powers of the metric perturbation hij. Next we turn to a discussion of the above 
results in different gauges. 



4 Gauge choices and corresponding transformations 



The previous discussion and summary focused exclusively on the comoving gauge choice for the 
metric, implicit in the definition of Eq. ()2.6p . Next we will consider some additional gauges. In 
this paper we will specifically refer to three choices for the metric: the comoving, synchronous and 
conformal Newtonian forms. The first two are closely related to each other, and were used to obtain 
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part of the results presented in our previous work [151 [T8] , which was summarized in the previous 
section. Note that in our previous work [TS] we did not include the effects of a stress field s, since 
it was not necessary for the discussion of density perturbations; new terms arising from such a field 
are included below. The third form of the metric is the primary focus of the present discussion; 
the results obtained later on in this paper will either be derived for this metric, or transformed to 
it by relying on results obtained previously in the other gauges. 



4.1 Comoving, synchronous and conformal Newtonian gauges 



The comoving metric has the form 



with background metric 



For the fluctuation one sets 



g^,v = diag [-\,c? 



hoi — hio — 



and decomposes the remaining hy as 

hij(k,t) = c? 



1 



ki kj 



(4.1) 

(4.2) 
(4.3) 

(4.4) 



so that Tr(hij) = a} h. Besides the scale factor a, the metric is therefore parametrized in terms of 
the two functions s and h. 

On the other hand, in the synchronous gauge one sets again g^^, = g^v + h^^ now with back- 
ground metric 

= diag (-1,1, 1,1) . (4.5) 
For the fluctuation one sets again /iQj = /lio = and 



hij(k,t) = a 



(4.6) 



so that now Tr{hij) = c? hgync- Here, besides the overall scale factor a, the metric is parametrized 
in terms of the two functions r] and hgync- From a comparison of the two gauges (comoving and 
synchronous) one has 



2r] 



1 

~3 
13 



(h + s) 



(4.7) 



and 

hsync + 6r/ = -s . (4.8) 

Finally the conformal Newtonian gauge is in turn described by two scalar potentials ip and (j)- 
In this case the line element is given by 

dr^ = -g^,ydx^ dx" = { (1 + 2^) dt^ -{1-2^) dxidx'^ . (4.9) 

Therefore for the metric itself one writes again g^u = g^u + with g^xv = diag (—1, 1, 1, 1) as 
for the synchronous case, and furthermore /iQi = /ijo = as before, and now 

/ioo = a^(-2V) (4.10) 

hij = a^ (-2^) 6ij . (4.11) 

A suitable set of gauge transformations then allows one to go from the synchronous, or comoving, 
to the conformal Newtonian gauge [19]. 



4.2 Tensor box in the comoving gauge 

To compute higher order contributions from the hij 's appearing in the comoving gauge metric, 
one needs to expand G{0) in the various metric perturbations, 



G{a) = Go 



^+i^((5k)'""-^5k-°'"(M).(^)"%..: 



(4.12) 



where the superscripts (0) and (1) refer to zeroth and first order in this expansion, respectively. 
To get the correction of 0{h, s) to the field equations, one therefore needs to consider the relevant 
term in the expansion of (1 + 5G{0)/Go) T^^, 



This last form allows us to use the results obtained previously for the FLRW case, namely 

I. ' t,, = (4.14) 

with here 

= \Pvac{'t) + Pvac{t)] Uy + gf„yPvac{t) (4.15) 
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to zeroth order in h, and 

/Omc(t) = p{t) Pvacit) = W^ac p{t) • (4.16) 

Go Go 

and Wvac = 1/3. Therefore, in light of the results of Ref. [Ej, the problem has been reduced to 
computing the more tractable expression 

-^^■o(^Hh,s).T;:^. (4.17) 

To make progress, we will assume a harmonic time dependence for both the perturbations h{t) = 
/lo e*"^* and s(t) = soe*'^*, and for the background quantities a{t) = aoe*^*, p{t) = p^e^^^, and 
6G{t) = 6Gq e*^*. From now on we shall consider both uj and T as slowly varying functions (indeed 
constants), with the time scale of variations for the perturbation much shorter than the time scale 
associated with all the background quantities. A more sophisticated treatment will be reserved for 
future work. Therefore we will take here uj^T or h/h^ a/ a, which is the same approximation 
that was used in obtaining the results of Ref. [18j . 

Let us now list, in sequence, the required matrix elements needed for the present calculation. 
For the tensor box tt matrix element {—-^ -^oj ■ O^'^'^^h, s) ■ T^"'^)oo one obtains 

1 11 5G(t) , , r , ,^,,2, , , 

ZU 6 Gq uj 

For the tensor box ti matrix element 5^07 " s) ■ T^°'^)Qi one obtains 

-ih^l^-^ pit) —{h-2s) + 0{k') . (4.19) 

For the tensor box ii matrix element, summed over i, {—-^ -^oj ■ □^^^(/i, s) • T^°''^)ii, one obtains 

3 (+ ^ w.ac ^ ^ Pit) ^h)+ Oik^) . (4.20) 
\ zv 3 Go J 

For the tensor box ii matrix element, not summed over z, (—27; 5^ " n^^^(/i, s) ■T'"°''^)ii, one obtains 



1 2 ^Git) , , 

+ 2;;'" IvT"*" 



Wvac + ^ (1-3t^ ) — S 

00 



3 a; 9 V k'^ 



+ Oik'^). (4.21) 



Finally for the tensor box ij matrix element, (— 217 " n^^^(/i, s) • T'""''^)ij, one obtains 



^^«i^pW + (4.22) 

fe^ ziy 6 Go UJ 



The above expressions are now inserted in the general effective field equations of Eq. (j2.4p , and will 
give rise to a set of effective field equations appropriate for this particular gauge, to first order in 
the field perturbation and with the effects of G'(n) included. 
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4.3 Field equations in the comoving, synchronous and conformal Newtonian 
gauges 



As a result of the previous manipulations one obtains in the comoving gauge with fields (/i, s) 
the following tt, ti, ii (or xx + yy + zz), and ij field equations 



3 a 



a 



Go 2v 



i(i + s)=8,G„|(-±)|l(.-2s)p + 0(t 



2v I 9 iu 



— ^ ih + s) -3-h-h = 247rGo — w^acpS + 24ttGo — w^ac ^,hp + 0{k^ 



3a^ 



a 











2v 



- -:t (n + s) - - -s - - s 



6 ' '2a 2 
As in Ref. |18j . we have found it convenient to here to set in the above expressions 

8\ r 
3/ w 



and 



6 LiJ A bJ 



(4.23) 
(4.24) 
(4.25) 
(4.26) 

(4.27) 
(4.28) 



In the field equations listed above the terms 0{k'^) arise because of terms 0{k^) in the expansion 
of the tensor box operator. 

The next step is to convert the left-hand sides of the above field equations, namely Eqs. (j4.23p . 
(j4.24p . (j4.25p and (|4.26p . which are all expressed in the comoving gauge {h,s), to the synchronous 
gauge with fields {hsync, f])- The result of this change of gauge is the sequential replacement 



{h + s) + -h 



k' 



1 d ■ 



3 a 



2 nV + 9 hgync 



2r] 



3^ 

1 e 



„ ^ (h + s) -3- h - h 
3 a 



k^ 

2 — 7] 



6 a 



2 {h + s) 



3 a 1 

s s 

2 a 2 



k^ 

T2^ + 



1 

1 1 



1 d ■ 



(/isync + 6 7/) + ^ - (As?/nc + 6 (4.29) 



a— 2 a2 

The next step involves one more transformation, this time from the synchronous {hgync-, v) to the 
desired conformal Newtonian ((/>, ip) gauge, 



-2k^r] 



''sync 



+ 3' 



+ 

a 



2r] 



+ - V 
a 
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0? 



a2 



aV ^ \ a a'^ J 3 

^(</.-V). (4.30) 



a2 



k^'n+ \ {hsync + 67?)+- (/lsy„c + 6 7? 



a? 



Equivalently, the above sequence of two transformations can be described by a single transforma- 
tion, from comoving (/i, s) to conformal Newtonian (0, tp) gauge, which is trivially obtained by 
combining the previous two. The final outcome of all these manipulations is to achieve a rewrite 
of the full set of four original field equations, given in Eqs. (j4.23p . (j4.24p . ()4.25p and ()4.26p . now 
with the left hand side given in the conformal Newtonian gauge and the right hand side left in the 
original comoving gauge. One obtains 

A;20 + 3^ (^ + ^^\=-AT,Goa' f 1 + ^"l p 6 - At^G^ a" ^-^ ^ h p + Oik") (4.31) 



a \ a J \ Go J Go 2i/ 

d \ . „ 5G f 1 \ 2 1 „ , _ „.,2n 



-+-7/; =4^Go— -- --{h-2s)p+ 0{k') (4.32) 
a J Go \ Zu J 9 luj 

+ 47rGo a — w^ac tt P 
Go Zv 

+ 0{k^) (4.33) 

A:2 _ ^) = + SttGo a' -^^sp 0{k^) . (4.34) 

Go Zu 

Note that we have, for convenience, multiplied out the first, third and fourth equations by a factor 
of a^. The last equation involves the quantity 

"=3G^2^-^- 

For the purpose of computing the gravitational slip function rj = Tp / (j) — 1 it will be useful here to 
record the following relationship between perturbations in the comoving and conformal Newtonian 
gauge. One has 

1 2 f ■■ -^0. 



Use has been made here of the following relationship between derivatives of an arbitrary function 
/ in the synchronous and comoving gauges 

jsync ^ ^jcom ^-^ ^g^ 
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and 

d d , ^ 

" ' sync ' com 

SO that 

f ^com \ 
jsync ^ ^2 / ^ jcom _|_ jcom \ ^ (4.40) 



5 Gravitational slip function 



The gravitational slip function is commonly defined as 

V-^. (5.1) 

In classical GR one has (p = ip so that ?7 = 0, which makes the quantity r] a useful parametrization 
for deviations from classical GR, whatever their origin might be. Using the ij field equation given 
in Eqs. (j4.3ip . (j4.32p . (j4.33p and (j4.34p . and the relationship between the conformal Newtonian 
fluctuation (j) and the comoving gauge fluctuations h and s, one finally obtains the rather simple 
result 

7] = —— = - levrGo — — - T P . (5.2) 
(p Gq Zu a s 



The last expression contains the quantity 



8\ iV 



(5.3) 



where ojs is the frequency associated with the s perturbation, and we have made use of — )• d/a. 
An equivalent form for the expression in Eq. ()5.2p is 

6G I S I s _ 1 8 /sdt _ 

7? = -IQttGq — — - — -p = -levrGo — — -^^p. (5.4) 

Go 2,v 6 lOJs s Gq zi' 6 s 

In the last expression we now can make use of the equation of motion for the perturbation s{t) to 
the order we are working, namely 

s + 3-s = . (5.5) 

a 

Let us look here first at the very simple limit of A ~ 0; the physically more relevant case of 
nonzero A will be discussed a bit later. Note that, in view of Eq. (12. 2p . this last limit corresponds 
therefore to a very large ^. Then for a perfect fluid with equation of state p = wp one has simply 
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a{t) = ao{t/tof/^^^+'^'^ and p{t) = l/[6TrGt'^{l + wf], and from Eqs. ([531) or (f5T71) one obtains for 
w = 



t\\ (t 



r? = 4.^Q (-j \^y^^]^0(f) (5.6) 



whereas for w 7^ one has 

^ = . (7) (5.7) 

3 w (1 — \^/ 

Another extreme, but nevertheless equally simple, case is a pure cosmological constant term (no 
matter of any type), which can be modeled by the choice w = —\. In this case t is related to the 
scale factor by 

a{t) 




exp<^ W-(t-to) ) ■ (5. 



Then, using the relation in Eq. (j2.2p . one obtains 

I = 1 + ln-, (5.9) 
where the quantity is therefore related to the time to ("today", oq = 1) and the scale ^ by 

^ = 1 + ln — . (5.10) 

Since numerically to is close to, but smaller than, ^, the scale factor will be close to, but slightly 
larger than, one. 

To actually come up with a definite number for r] in more realistic cases, one needs (apart 
from including the effects of A 7^ 0, which is done below) a value for the coefficient q appearing 
in Eq. (|2.1ip for G{t), which in turn is related to the original expression for the running Newton's 
constant G(n) in Eq. ()2.ip . This issue will be discussed in some detail later, but here let us say the 
following. In Ref. |15j it was estimated that the values of q in Eq. (|2.1ip and cq in Eq. (|2.1ip are 
of the same order of magnitude, ct ~ 0.62co. The most difficult part has been therefore a reliable 
estimate of cq, which is obtained from a lattice computation of invariant curvature correlations 
at fixed geodesic distance [20] . and which, after reexamination of various systematic uncertainties, 
leads to the recent estimate used in|18] of cq ~ 33.3. That would give q ~ 20.6 which, as we will 
see later, is still very large. Nevertheless it is expected that cq (or q) enter all calculations with 
G(n) with the same magnitude and sign. 

Let us now go back to the more physical case of A 7^ 0. The relevant expression for 7]{t) is 
Eq. ()5.4p . where we use the equation for s(t), Eq. ()5.5p . to eliminate the latter. It is also convenient 
at this stage to change variables from t to a{t), and use the equivalent equation for s(a), namely 
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where the prime denotes differentiation with respect to the scale factor a. In the above equation 
one can use, for nonrelativistic matter with equation of state such that w = and to the order 
needed here, the first Friedmann equation 



We have also made use of the unperturbed result for the background matter density valid for w = ^ 
(which follows from energy conservation), namely 

P = Po\- (5-13) 

Note that the above expression for p is valid to zeroth order in 5G, which is entirely adequate when 
substituted into /?(a), since the rest there is already first order in 5G. This finally gives an explicit 
solution for s(a) 

s{a) oc -^T^ Vl + a3^, (5.14) 

with parameter 9 = A/SvrGo po- The above solution for s(a) can then be substituted directly 
in Eq. (j5.4p . provided one changes variables from t to a(t), and in the process uses the following 
identities 



sit) At = sia) ——- da , (5.15) 
J aid [a) 

as well as 

s = aH(a)—, (5.16) 
oa 

with H{a) given a few lines above. 

The resulting expression, which still involves an integral over the scale factor a(t), can now be 
readily evaluated, and leads eventually to a rather simple expression for r]. The general result for 
nonrelativistic matter (w = 0) but A / is 



16 6G{a) 

ma) = IT- —p^ 



a 



(5.17) 



This is the main result of the paper. The integration constant has been fixed following the 
requirement that the scale factor a — ?• for t — )• ^ [see Eqs. ()2.ip . (12. lip and (j3.8p for the definitions 
of In other words, by switching to the variable a{t) instead of t, the quantity ^ has been traded 
for a^. In the next section we will show that in practice the quantity is generally expected to be 
slightly larger than the scale factor "today", i.e. for t = to- As a result the correction in Eq. (|5.17p 
is expected to be negative today. 
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The next section will be devoted to establishing the general relationship between t and a{t), for 
nonvanishing cosmological constant A, so that a quantitative estimate for the slip function r/ can 
be obtained from Eq. (j5.17p in a realistic cosmological context. Specifically we will be interested 
in the value of rj for a current matter fraction Q ~ 0.25, as suggested by current astrophysical 
measurements. 



5.1 Relating the scale factor a to t, and vice versa 



Let us now come back to the general problem of estimating 77(a), using the expression given 
in Eq. (|5.17p . for A 7^ and a nonrelativistic fluid with w = 0. To predict the correct value for 
the slip function 77(a) one needs the quantity 6G{a), which is obtained from the FLRW version of 
G(n), namely G{t) in Eq. (|2.1ip . via the replacement, in this last quantity, of t — t- t{a). The last 
step requires therefore that the correct relationship between t and a{t) be established, for any value 
of A. In the following we will first relate t to a{t), and vice versa, to zeroth order in the quantum 
correction 6G [we will call them a^^\t) and t^'^\a)], and then compute the first order correction in 
6G to the above quantities [we will call those a^^\t) and t^^\a)]. 

Let us look first at the zeroth order result. The field equations and the energy conservation 
equation for a^'^\t), without a 5G correction, but with the A term, were already given in Eq. (|2.17p . 

3^ = 8vrGop-(°)(t) + A 

a(0)2(t) aW(t) ^ mu^ . , ^ 

-^ + 2-^ = -8-Go^.-(°)(t) + A (5.18) 

for a spatially flat universe (A; = 0), and 

/)(t) + 3(l + u;)^p-W(t) = 0. (5.19) 

From these one can obtain a^^\t) and then p^^\t). As a result the scale factor is found to be related 
to time by 



t(°)(a) = ^ (5.20) 



2 Arcsinh 



,3/2 



/3A 

where we have deflned the parameter 

A 1-n 
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with pq the current (t = to) matter density, and 17 the current matter fraction. Note that in practice 
we win be interested in a matter fraction which today is around 0.25, giving 6 ~ 3.0, a number 
which is of course quite far from the zero cosmological constant case of ^ = 0. 

One can express the time today (to) in terms of cosmological constant A, and therefore in terms 
of 0, as follows 



t. 



(0) 



2 Arcsinh (\/^) 



(5.22) 



with the normalization for t^^\a) such that t^^\a = 0) = and t^^\a = 1) = to "today". So here 
we follow the customary choice of having the scale factor equal to one " today" . Then one has 

Arcsinh 



tW(a) 

,(0) 



t^^' Arcsinh(\/^) 

When expanded out in 9, the above result leads to some perhaps more recognizable terms, 

1 



(5.23) 



t(0)(a) 

— a2 



(0) 



1 



+ 



-(- 

560 V 



f-l + a , - , 
6 V / 360 

Conversely, one has for the scale factor as a function of the time 

/Sinh^ ^ 



17- lOa^ + 270*^1 + 



(5.24) 



aW(t) 



(5.25) 



which, when expanded out in A or t, gives the more recognizable result 



a(°)(t) 



Similarly for the pressure one obtains 



3At2 / At2 A^t^ ' 



(5.26) 



p(0)(t) 



A Csch^ 



3A 



SvrGo 



which when expanded out in A or t gives the more familiar result 



p-W(t) 



1 



6vrGotMl + ^ + 4? + ^ + 



(5.27) 



(5.28) 



To be more specific, let us set ^ = 3, which corresponds to a matter fraction today of ~ 0.25. In 
addition, we will now make use of Eq. (|2.2p and set A — t- 3/^^. One then obtains 



.(0) 



3) = 0.878 ^ , 



(5.29) 



which shows that to and ^ are rather close to each other (apparently a numerical coincidence). 
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Then, from the expression for G{t) in Eq. (I2.1ip . 



1 



*°'*'-c.m-. (5.30) 



Go V?, 

one can obtain G{a) in all generality, by the replacement t — >• t{a) according to the result of 
Eqs. (|5.20p or (j5.23p . For the special case of pure nonrelativistic matter with equation of state 
w = and A = one obtains, using Eq. ()5.24p . 

'-^ - . i-Y . (5.31) 

Go \a!:J 

with exponent 

= (5.32) 

The latter is largely the expression used earlier in the matter density perturbation treatment of 
our earlier work of Ref. [18j. 

More generally one can define as the value for the scale factor a which corresponds to the 
scale ^, 



af = l^] Sinht 



1.655 (l)' , (5.33) 



so that in general / oq, where ao = 1 is the scale factor "today". Then for the observationally 
favored case ^ ~ 3 one obtains 

af\e = 3) = 1.U8 , (5.34) 

which clearly implies a^^'' > oq = 1. 3 The above expressions will be used in the next section to 
obtain a quantitative estimate for the slip function 77(a), evaluated at today's time t = Iq. 

The discussion above dealt with the case of 5G = 0. Let us now consider briefly the corrections 
to a{t) and, conversely, t{a) that come about when the running of G is included, in other words 
when a constant G is replaced by G{t) or G{a) in the effective field equations. In Eq. (|2.18p the 
Friedman equations were given in the presence of a running G, namely 



2m+2-7TT = -8ttGo [w + w^ac^r^] pit) + >^, (5-35) 
a^(t) a[t) \ Go J 



together with the energy conservation equation 



{l + w) + {l + w, ' 



a{t) r ' ' - ' -^-^ Go 



p-W + ^p-W+(l + ^)Mt) = 0. (5.36) 



Let us give here a few more observational numbers for present and future reference. From the present age of 
the Universe to « 13.75 Gj/rs ~ 4216 Afpc, whereas from the observed value of A (mostly extracted from distant 
supernovae surveys) one has following Eq. 1)2. 2[) ^ ~ 4890 Mpc, which then gives to/£, — 0.862 = 1/1.160. This last 
ratio is similar to the number we used in Eq. (|5.29p . by setting there Q — 0.25 exactly. 
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To solve these equations to first order in 5G we set 

a{t) = a^'^\t)[l + cta^^\t)\ (5.37) 

/^(t) = /^W(t)[l + QP«(t)] (5.38) 

where a*-''-*(t) and p^^\t) here represent the solutions obtained previously for 5G = 0. One then 



finds for the correction to the matter density 



p-(^)(i) 



and to lowest nontrivial order in t and for w^ac 

p-W(t): 



(1 + ^) 

= 1/3 



3 + 5z^ fr ^ 



'2,\t Coth 



/3A 



3(1 + Z.) ve 

For the correction to the scale factor one finds 



a«(t) 



A 



{l + v) Jo _i + Cosh 
and to lowest nontrivial order in t for Wyac — 

1/3, 



/3Xt' 



dt' 



aW(t) 



t ^ 



(5.39) 



(5.40) 



(5.41) 



, , ^ , (5.42) 

9(l + z.) \CJ 

After having obtained the relevant formulas for a{t) and t{a) in the general case, i.e. for nonzero 
A, we can return to the problem of evaluating the slip function rj. 



5.2 Quantitative estimate of the slip function 77(2;) 



The general expression for the gravitational slip function r]{a) was given earlier in Eq. (I5.17P 
for w = and A / 0, 

(5.43) 



16 5G{a) 
via) = — log 



3i^ Go 

To obtain 5G(a) we now use, from Eq. (|2.11[) . 

5G{t) 



a 



Go 



(5.44) 



and substitute in the above expression for 6G{t) the correct relationship between t and o, namely 
t{a) from Eq. (|5.20p . which among other things contains the constant defined in Eq. (|5.33p . 



1\ 3 2 
-1 Sinh3 



(5.45) 
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It will be convenient, at this stage, to also make use of the relationship in Eq. (12. 2p . namely 

A ^ ^ . (5.46) 

The last step left is to make contact with observationally accessible quantities, by expanding in the 
redshift z, related in the usual way to the scale factor a by a = 1/(1 + z). Then for i' = 1/3 and 
6 = 3 (matter fraction 0, = 0.25) one finally obtains for the gravitational slip function 

r]iz) = -1.491 ct -6.418 C£Z + 30.074c£z2 + ••■ (5.47) 

To obtain an actual number for r]{z = 0) one needs to address two more issues. They are (i) to 
provide a bound on the theoretical uncertainties in the above expression, and (ii) to give an estimate 
for the coefficient q, which is traced back to Eq. (|2.1ip and therefore to the original expression for 
G(n) in Eq. (12. ip . The latter contains the coefficient cq, but in Ref. [15j the estimate was given 
Ct = 0.450 Co for the tensor box operator; thus q and cq can safely be assumed to have the same 
sign, and comparable magnitudes. 

To estimate the level of uncertainty in the magnitude of the correction coefficient in Eq. (j5.47p 
we will consider here an infrared regulated version of G(n), where an infrared cutoff is supplied 
so that in Fourier space k > and the spurious infrared divergence at small k is removed. 
This is quite analogous to an infrared regularization used very successfully in phenomenological 
applications to QCD heavy quark bound states [2H I22j. and which has recently found some limited 
justification in the framework of infrared renormalons [23]. As shown already in the first cited 
reference, it works much better than expected; here a similar prescription will be used just as a 
means to provide some estimate on the theoretical uncertainty in the result of Eq. ()5.47p . Therefore, 
instead of the G{0) in Eq. (|2.ip . which in momentum space corresponds to 



G(r) ~ G| 

we will consider a corresponding infrared regulated version 



1 \ 1/21^ 



Gi^) ^ Go 



k^' + C 



2 \ 1/2^ 



1 + CO I TT—J^ + 



(5.48) 



(5.49) 



Of course the small distance, k ^ S^^^ ov r <^ ^, behavior is unchanged, whereas for large distances 
r ^ ^ the gravitational coupling no longer exhibits the spurious infrared divergence; instead it 
approaches a finite value Gqo — (1 + cq + . . .) Go- Now, in momentum space the infrared regulated 
6G{k) reads 

j^r\ (-0, 



Go \ + m 
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with m = 1/^, and in position space the corresponding form is 



6G(a) / 1 y/^'' 

- CO ^^^-r • (5.51) 



Following the results of Ref. [15] , if the above differential operator acts on functions of t only, then 
one obtains for 5G{t) 



6G{t) 

-p=r- = Co 



/ 



1 
2 1^ 



1 



(5.52) 



with again q/co ~ 0.62 [15] . Note that the expression in Eq. (j5.52p could also have been obtained 
directly from Eq. (|2.11|) . by a direct regularization. 

One can then repeat the whole calculation for ry(a) with the regulated version of 5G{t) given in 
Eq. ([532]) . The result is 

r^[z) = -0.766 ct - 4.109 q z + 12.188 ct + ■ ■ ■ . (5.53) 

It seems that the effect of the infrared regularization has been to reduce the magnitude of the 
effect (at z = 0) by about a factor of 2. It is encouraging that, at this stage of the calculation, the 
negative trend in r]{z) due to the running of G appears unchanged. Furthermore, in all cases we 
have looked so far, the value r]{z = 0) is found to be negative. 



5.3 Slip function 77(2;) for stress perturbation s = 

In Ref. [18j a preliminary estimate of the magnitude of the slip function rj was given. The 
calculation there neglected the stress field s in Eq. (j4.4p and only included the metric perturbation 
h in the comoving gauge. The main reason was that nonrelativistic matter density perturbations, 
and therefore the growth exponents, are unaffected by the stress field contribution. We will show 
here that in this case one still obtains a nonvanishing r], whose value we will discuss below. The 
results will be useful, since now a direct comparison can be done with the full answer (including 
the stress field) for r]{z) given in the previous section. 

In the absence of stress (s = 0) and finite k, the tt and xx -\- yy + zz field equations read 

-2^^- SttGo^ ^ p S (-j^) = 8vrGo f 1 + ^) (5.54) 
a"^ 2i/ Go \ l + wj \ Go J 

2%. (-0 - (A) + 247rGo^ w^ac ( —r^) = - 247rGo (w + w^ac ^) pS ■ (5.55) 



2v Gq \ \ ^ w ) V Gq 
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In both equations we have made use of zeroth order (in 6G/Gq) energy conservation, which leads 
to h = —jj^:^ S, where 6 is the matter fraction. One can then take the ratio of the two equations 
given above, and obtain again an expression for the shp function rj = (^ — 0) /(j). For w = 
(nonrelativistic matter), after expanding in 6G/G0, one finds the rather simple result 

Here the quantity is the same as in Eq. (j4.28p . and depends on the choices detailed below. In 
the following we will continue to use Wvac = 1/3 [see Eqs. ()2.15p and (12.161) ] [El US], which is the 
correct value associated with G(n) in the FLRW background metric. 
In Ref. [H] we used the scalar box value Ch = 1/2, which then gives 

In this last case it is then easy to recompute the slip function in terms of the redshift, just as was 
done in the previous section, and one finds, under the same conditions as before [v = 1/3, = 3, 
and to/C as given in Eq. (|5.29p ] the following result 

ry ~ - 0.338 Q + 0{z) . (5.58) 

For the infrared regulated version of 6G/Gq given in Eq. (I5.52P one obtains instead the slightly 
smaller value 

r] ~ - 0.174 ct + 0{z) . (5.59) 

For the tensor box case (also discussed extensively in fTS], where it was shown that this is in fact 
the correct way of doing the calculation) one finds a significantly larger value — 7.927, so that 
in this case the slip function 77 becomes 



, 7.927\ 5G ( 7.927\ ( t\- , 

Also in this case one can recompute the slip function in terms of the redshift, and one finds, under 
the same conditions as before, 

r/~ - 15.42 Q + 0{z) . (5.61) 
For the infrared regulated 6G/Gq given in Eq. (15.521) one finds instead 

r] ~ - 7.919 Q + 0{z) , (5.62) 

which is again about a factor of 2 smaller than the unregulated value. 
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We conclude from the above exercise of calculating ij with vanishing stress field s = three 
things. The first is that using the scalar box result on the trace of the energy-momentum tensor 
(which ultimately is not an entirely correct, or at least an incomplete, procedure, given the tensor 
nature of the matter energy- momentum tensor) underestimates the effects of (?(□) on the slip 
function r]{z = 0) by a factor that can be as large as an order of magnitude. 

The second lesson is that the stress field (s) contribution is indeed important, since it reduces 
the size of the quantum correction significantly [Eqs. ()5.47p and ()5.53p ]. compared to the s = 
result [Eqs. (I5.60p and (I5.61j) ]. again by almost an order of magnitude, which would imply some 
degree of cancellation between the s and h contributions. 

The third observation is that in all cases we have looked at so far the quantum correction to 
the slip function is negative at z = 0. 

6 Conclusions 

In the previous sections we computed corrections to the gravitational slip function r] = — l 
arising from the renormalization-group motivated running (?(□), as given in Eq. (j2.ip . The relevant 
result was presented in Eqs. (j5.47p and (j5.53p . the first expression representing the answer for an 
unregulated G'(n), and the second answer found for an infrared regulated version of the same. It 
should be noted that, so far, in the treatment of metric and matter perturbations we have considered 
only the k — ^ limit [see Eq. ()3.2p ]. Let us focus here for definiteness on the first of the two results 
[Eq. dOTIl ]. which is 

r}{z) ~ -1.491 Q + 0{z) (6.1) 

at z ~ 0. We now come to the last issue, namely an estimate for the magnitude of the constant 
ct- As already discussed previously in Sec. 15. 2^ to get an actual number for rj{z = 0) one needs 
a number for q, whose appearance is traced back to Eq. (|2.1ip . and therefore to the original 
expression for G(n) in Eq. (j2.ip . with q ~ 0.450 x cq for the relevant tensor box operator [15j . 

The value of the constant cq has to be extracted from a nonperturbative lattice computation 
of invariant curvature correlations at fixed geodesic distance [20]; it relates the physical correlation 
length ^ to the bare lattice coupling G, and is therefore a genuinely nonperturbative amplitude. 
After a reexamination of various systematic uncertainties, these lead to the recent estimate used 
in [18] of Co ~ 33.3. That would give for the amplitude q ~ 20.6 which still seems rather large. 
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Nevertheless, based on experience with other field-theoretic models which also exhibit nontrivial 
fixed points such as the nonlinear sigma model, as well as QCD and non-Abelian gauge theories, 
one would expect this amplitude to be of order unity; very small or very large numbers would seem 
rather atypical and un-natural. 

As far as astrophysical observations are concerned, current estimates for r]{z = 0) obtained from 
CMB measurements give values around 0.09ib0.7 |301l31j. which would then imply an observational 
bound ct < 0.3. 

Indeed a similar problem of magnitudes for the theoretical amplitudes was found in our recent 
calculation of matter density perturbations with G'(n), where again the corrections seemed rather 
large [18] in view of the above quoted value of q. Let us briefly summarize those results here. 
Specifically, in Ref. [18] a value for the density perturbation growth index 7 was obtained in the 
presence of G{0). The quantity 7 is in general obtained from the growth index /(a) [2l] 



_ din 6{a) 

JW = — , (6.2) 

dm a 

where 5{a) is the matter density contrast. One is mainly interested in the neighborhood of the 
present era, a{t) ~ oq — 1, which leads to the definition of the growth index parameter 7 via 

In/ 



7 = 



(6.3) 

a=aQ 



Inn 

The latter has been the subject of increasingly accurate cosmological observations, for some recent 
references see [251 [Ml 127] . |f 

On the theoretical side, for the tensor box one finds [18], for a matter fraction = 0.25, 



7 = 0.556 - 106.4 Q + 0(c?) . (6.4) 

where the first contribution is the classical GR value from the relativistic treatment of matter 
density perturbations [24J. The result presented above is in fact a slight improvement over the 
answer quoted in our earlier work [IS], since now the improved relationship between t and a given 
in Eq. ()5.20p has been used, which reduces the magnitude of the correction proportional to q. 
Nevertheless, it should be emphasized that the above result has been obtained in the k — )■ limit 
of the perturbation Fourier modes in Eq. (|3.2p . 

Recent observational bounds on x-ray studies of large galactic clusters at distance scales of up 
to about 1.4 to 8.5AIpc (comoving radii of ~ 8.5Mpc and viral radii of ~ lAMpc) [26j favor values 
for 7 = 0.50 lb 0.08, and more recently 7 = 0.55 + 0.13 — 0.10 [27j. This would then constrain the 



^ For a recent detailed review on the many tests of general relativity on astrophysical scales, and a much more 
complete set of references, see for example [281 129| . 
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amplitude q in Eq. (16. 4p at that scale to ct ^ 5 x 10"'^. The latter bound from density perturbations 
seems a much more stringent bound than the one coming from the observed slip function. Indeed 
with the bound on q coming from the observed density perturbation exponents one would conclude 
that, according to Eq. (16. ip . the correction to the slip function at z ~ must indeed be very small, 
r] ~ 0(10^'^), which is a few orders of magnitude below the observational limit quoted above, 
r] ~ 0.09 ± 0.7. 

It is of course possible that the galactic clusters in question are not large enough yet to see the 
quantum effect of G(n), since after all the relevant scale in Eq. (12. ip is related to A and is supposed 
therefore to be very large, ^ ~ 4890Mpc. |^ But most likely the theoretical uncertainties in the 
value of Ct have also been underestimated in |20J, and new, high precision lattice calculation will 
be required to significantly reduce the systematic errors. 

Nevertheless it seems clear that the non-perturbative coefficient cq (or q) enters all calculations 
involving G(n) with the same magnitude and sign. This is simply a consequence of cq being part 
of the renormalization group G{0) which enters the covariant effective field equations of Eq. (j2.4p . 
Consequently, one should be able to relate one set of physical results to another, such as the value 
of the slip function ri{z = 0) in Eq. (|5.47p to the corrections to the density perturbation growth 
exponent 7 computed in [18j, and given here in Eq. 1)6.41) . Then the amplitude q can be made 
to conveniently drop out when computing the ratio of G{0) corrections to two different physical 
processes. The resulting predictions are then entirely independent of the theoretical uncertainty 
in the amplitude cq, and remain sensitive only to the uncertainties in the two other quantum 
parameters ^ and z^, which are expected to be significantly smaller. One then obtains for the ratio 
of the corrections to the growth exponent 7 to the slip function r]{z = 0) at z ~ 

(57 -106.4Q , 

~ ~ +71.4 (6.5 

6rj -1.491 Q ^ ^ 

for the infrared unimproved case. One conclusion that one can draw from the numerical value 

of the above ratio is that it might be significantly harder to see the G{0) correction in the slip 

function than in the matter density growth exponent, by almost 2 orders of magnitude in relative 

magnitude. Hopefully increasingly accurate astrophysical measurements of the latter will be done 



® One might perhaps think that the running of G envisioned here might lead to small observable consequences on 
much shorter, galactic length scales. That this is not the case can be seen, for example, from the following argument. 
For a typical galaxy one has a size ~ 30 kpc, giving for the quantum correction the estimate, from the potential 
obtained in the static isotropic metric solution with G{0) [16] which gives 5G{r) ~ (r/O^/", (30fcpc/4890xl0^ kpcf ~ 
2.31 X 10~^^ which is tiny given the large size of ^. It is therefore unlikely that such a correction will be detectable at 
these length scales, or that it could account for large anomalies in the galactic rotation curves. The above argument 
also implies a certain sensitivity of the results to the value of the scale ^ ; thus an increase in ^ by a factor of two tends 
to reduce the effects of G{D) by roughly 2^ = 8, as can be seen already from Eq. (|2.ip with v — 1/3. More specifically, 
the amplitude of the quantum correction is proportional, in the non-infrared improved case, to the combination co/^^. 
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in the not too distant future. Of particular interest would be any trend in the growth exponents 
as a function of the maximum galactic cluster size. 
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Appendix 

A Scalar box in the comoving gauge 



In this section we will give a short sample calculation of the effects of the covariant d'Alembertian 
operator □ = g^'^ 'V^VnU acting on a coordinate scalar, such as the trace of the energy-momentum 
tensor. The calculation presented below will show that the result is unchanged when the stress 
contribution s is included in the metric for the comoving gauge. Specifically here we will be 
interested in the correction of order hij that arises when the operator in Eq. (14.121) acts on the 
scalar T^'^ = —p- Thus, for example, it will give the correction 0{h, s) to 6pyac, namely the second 



term in the expression 



JG(aW) 6G{a){h,s) 

OPvac{t> = 7; dp{t) H p{t) 



(A.l) 



Go Go 

with the first term being simply given in the FLRW background by 6G{t)/GQ ■ dp{t). Here the 
0{h, s) correction is given explicitly by the expression 



SG{n){h,s) 



1 CO 1 



1 \l/2v 



Go 2u i^l^ 0(0) ' ' ' \ui^\ 

Now the covariant d'Alembertian □ acting on general scalar functions S{x) simplifies to 

In the absence of hij fluctuations this gives for the metric in the comoving gauge 

a(o)s{x) = \v^S-3^S-S . 



(A.2) 



(A.3) 



(A.4) 



To first order in the field fluctuation hij of the comoving gauge one computes 



S 



2 



1 



3^(' + ^) + il^ 



+ OAS 



2 hx ky 
a? fc2 



where we have set as usual h{x) = h{t) e^^'^. But, for a function of time only, one obtains 



1 



h{t)S{t) 



(A.5) 



(A.6) 



Thus to first order in the fluctuations one has 
1 



□ (0) 



• aW{h) ■ {6Gp) 



-df -3^dt 

t a 
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\hl^'^5G-5G^p (A.7) 



and there is no change from the result quoted in There we set s = 0, since we were only 

interested in cosmological density perturbations 5, which couple only to the trace part of the 
gravitational field fluctuations hij. 



34 



References 



[1] T. Damour, in Review of Particle Physics, J. Phys. G 33, 1 (2006); update in 



http://pdg.lbl.gov/2009/reviews/rpp2009-rev-gravity-tests.pdf (Nov. 2009). 



[2] T. Damour and K. Nordtvedt, Phys. Rev. Lett. 70, 2217 (1993); 
T. Damour and G. Esposito-Farese, Phys. Rev. D 53, 5541 (1996). 



[3] G. Veneziano, "String cosmology: The pre-big bang scenario," hep-th/0002094 



M. Gasperini and G. Veneziano, |hep-th/0703055 



[4] K. G. Wilson, Phys. Rev. Lett. 28, 548 (1972); Phys. Rev. D 7, 2911 (1973); see also: 

G. Parisi, in the proceedings of the 1976 Cargese Summer Institute, vol. 26 (Plenum Press, 
New York 1977); 

E. Brezin and J. Zinn-Justin, Phys. Rev. Lett. 36, 691 (1976); Phys. Rev. D14, 2615 (1976). 

[5] S. Weinberg, in 'General Relativity - An Einstein Centenary Survey', edited by S.W. Hawking 
and W. Israel, (Cambridge University Press, 1979). 

[6] J. Zinn-Justin, Quantum Field Theory and Critical Phenomena, Oxford University Press (New 
York, fourth edition, 2002). 

[7] H. W. Hamber, Quantum Gravitation, (Springer Publishing, Berlin and New York, 2009), and 
references therein. 

[8] H. W. Hamber and R. M. Williams, Nucl. Phys. B 248, 392 (1984); ibid. 260, 747 (1985); 
ibid. 269, 712 (1986); Phys. Lett. B 157, 368 (1985); Nucl. Phys. B435 361 (1995). 

[9] H. W. Hamber, Phys. Rev. D 45, 507 (1992); Nucl. Phys. B400, 347 (1993); Phys. Rev. D 
61, 124008 (2000). 

[10] H. Kawai and M. Ninomiya, Nucl. Phys. B336, 115 (1990); 

H. Kawai, Y. Kitazawa and M. Ninomiya, Nucl. Phys. B393, 280 (1993) and B404 684 (1993); 
Y. Kitazawa and M. Ninomiya, Phys. Rev. D55, 2076 (1997); 

T. Aida and Y. Kitazawa, Nucl. Phys. B491, 427 (1997). 

[11] M. Renter, Phys. Rev. D 57, 971 (1998); 

M. Renter, and H. Weyer, Gen. Relativ. Gravit. 41, 983 (2009); 

E. Manrique, M. Renter and F. Saueressig, hep-th/1006.0099, and references therein. 



35 



[12] G. A. Vilkovisky, in Quantum Theory of Gravity, edited by S. Christensen (Hilger, Bristol, 
1984); Nucl. Phys. B234, 125 (1984). 

[13] T. R. Taylor and G. Veneziano, Nucl. Phys. B345, 210 (1990); Phys. Lett. B228, 311 (1989). 

[14] H. W. Hamber and R. M. Wilhams, Phys. Rev. D 73 044031 (2006); ibid. D 76 084008 (2007); 
ibid. D 81 084048 (2010). 

[15] H. W. Hamber and R. M. Williams, Phys. Rev. D 72, 044026-1-16 (2005). 

[16] H. W. Hamber and R. M. Wilhams, Phys. Lett. B643, 228 (2006); Phys. Rev. D 75, 084014 
(2007). 

[17] D. Lopez Nacir and F. D. Mazzitelh, Phys. Rev. D 75, 024003 (2007). 

[18] H. Hamber and R. Toriumi, Phys. Rev. D 82, 043518 (2010). 

[19] C. -P. Ma and E. Bertschinger, Astrophysical Journal, v. 455, p. 7 (1995). 

[20] H. W. Hamber, Phys. Rev. D50 3932 (1994). 

[21] J. L. Richardson, Phys. Lett. B 82, 272 (1979). 

[22] E. Eichten and F. Feinberg, Phys. Rev. D 23, 2724 (1981); W. Buchmuller and S. H. H. Tye, 
Phys. Rev. D 24, 132 (1981); U. Ellwanger, M. Hirsch and A. Weber, Eur. Phys. J. C 1, 563 
(1998). 

[23] M. Beneke, Phys. Rept. 317, 1 (1999); M. Beneke and V. M. Braun, Nucl. Phys. B 454, 253 
(1995). 

[24] P. J. E. Peebles, Principles of Physical Cosmology, (Princeton University Press, 1993). 

[25] F. Schmidt, A. Vikhhnin and W. Hu, Phys. Rev. D 80, 083505 (2009). 

[26] A. Vikhlinin et al. larXiv:0903.5320l [astro-ph.CO]: la?Xiv:0903.2297l [astro-ph.CO]. 

[27] D. Rapetti, S. W. Allen, A. Mantz and H. Ebehng, larXiv:0909.3098l larXiv:0909.3099t 
larXiv:0911.1787 [astro-ph.CO]. 

[28] J. P. Uzan, Rev. Mod. Phys. 75, 403 (2003). 

[29] J. P. Uzan. .arXiv:0908.2243l [astro-ph.CO]. 



36 



[30] L. Amendola, M. Kunz and D. Sapone, JCAP 0804, 013 (2008). 
[31] S. F. Daniel et al, Phys. Rev. D 80, 023532 (2009). 



37 



